function [fig_h, edges_h] = DisplayEdgesByFrame(colony, fig_h, edges_h)

if nargin < 2, fig_h = figure; end

figure(fig_h);
hold on;
set(gca, 'xlimmode','manual',...
'ylimmode','manual',...
'zlimmode','manual',...
'climmode','manual',...
'alimmode','manual');

bounds = colony.bounds;
EBF = colony.edgesByFrame;
numImage = length(bounds);

for i=1:numImage
    if ~isempty(EBF(i).edges)
        edges_h(i) = plot(bounds(i)*ones(1, length(EBF(i).edges)) + (2*colony.margin(1)+1)/2, ...
            EBF(i).edges, 'r+');
    end
end

end